Phase-controlled superconducting heat-flux quantum modulator 
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We theoretically put forward the concept of a phase-controlled superconducting heat-flux quantum 
modulator. Its operation relies on phase-dependent heat current predicted to occur in temperature- 
biased Josephson tunnel junctions. The device behavior is investigated as a function of temperature 
bias across the junctions, bath temperature, and junctions asymmetry as well. In a realistic Al-based 
setup the structure could provide temperature modulation amplitudes up to ~ 50 mK with flux-to- 
temperature transfer coefficients exceeding ~ 125 mK/$o below 1 K, and temperature modulation 
frequency of the order of a few MHz. The proposed structure appears as a promising building-block 
for the implementation of novel-concept caloritronic devices operating at cryogenic temperatures. 



The investigation of thermal transport in solid-state 
nanodevices, usually referred to as caloritronics, is a 
research field which is nowadays garnering an increas- 
ing attention [TJ [5] since the impressive advances in 
nanoscience. The latter are suggesting that caloritron- 
ics is a central issue at the nanoscale, where heat plays a 
key role in determining the physical properties of the sys- 
tem. In this context, manipulation and control of heat 
currents is of great relevance for a number of applica- 
tions, and conceiving devices allowing to tune the amount 
of heat transferred through nanocircuits would represent 
an important breakthrough [51110]. Toward this direc- 
tion the prototype of a thermal modulator in the form 
of a heat transistor has been recently reported [HI [12] 
with the demonstration of electrostatic control of heat 
currents flowing through a superconducting circuit. 

Here we theoretically put forward the concept of 
a phase-controlled superconducting heat-flux quantum 
modulator. This device exploits the phase-dependent 
thermal current [13] predicted to occur in temperature- 
biased Josephson tunnel junctions to realize heat inter- 
ference. With a realistic design, temperature modula- 
tion amplitudes up to ~ 50 mK with flux-to-temperature 
transfer functions as large as ~ 125 mK/$o could be 
achieved below 1 K, and temperature modulation fre- 
quencies of the order of a few MHz. 

The structure we envision is sketched in Fig. [T] It con- 
sists of a double-tunnel-junction superconducting quan- 
tum interference device (i.e., a DC SQUID) composed of 
two superconductors Si and S2 in thermal equilibrium 
kept at different temperatures Ti and T 2 , respectively. 
Ra(b) denotes the normal-state resistance of Josephson 
junction J (m, while ip a and (ft, represent the macro- 
scopic quantum phase differences over the junctions. By 
neglecting the geometric inductance of the ring it fol- 
lows that if a + (fib + 27r<I>/$o = 2kir where <E> is the ap- 
plied magnetic flux through the loop, k is an integer and 
$0 = 2.067 x 10~ 15 Wb is the flux quantum. For definite- 
ness, we assume that Ti > T 2 so that the SQUID is only 
biased with a temperature drop across the junctions, but 
the voltage across them vanishes. In these conditions the 
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FIG. 1. (Color online) Scheme of the proposed device. Two 
superconductors Si and S2 kept at temperature Ti and Tb 
(with Ti > T2), respectively, are tunnel coupled so to imple- 
ment a DC SQUID. f a Q>) is the quantum phase difference over 
junction J a (j,) having normal-state resistance R a (b), is the 
applied magnetic flux threading the SQUID loop, and Qtot(<&) 
is the total heat current flowing through the structure. 

total heat current (Qtot) flowing from Si to S2 becomes 
stationary [13], and can be written as [T514T7] 

Qtot = Q qp (Ti,T 2 ) - Qint(T 1 ,T 2 ,(p a ,ipb)- (1) 

In Eq. 0, 

Q qp = Q a qp {T 1 ,T 2 )+Q\ p {T u T 2 ) (2) 

is the usual total thermal flux carried by quasiparticles 
through both junctions J a and J;,. On the other side, 

Qint = Qt n t(TuT 2 )cosip a + Q\ nt {T u T 2 )cos^, (3) 

describes the sum of the heat currents due to interfer- 
ence between quasiparticles and Cooper pairs condensate 
flowing through each junction of the SQUID. We empha- 
size that Qi n t originates from the Josephson effect and is 
peculiar to weakly-coupled superconductors [13TfT7] . Ac- 
cording to Eqs. ([!]) and (3), Q to t therefore consists of 
a phase-dependent function (Qi n t) superimposed on top 
of a phase-independent component. In writing Eq. ([I]) 
we have explicitly omitted the term related to heat cur- 
rent carried by the Cooper pairs condensate, as this term 
produces no dissipation in the system [131 HI HI] • 
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Aj(Tj) 2 ] is the normalized BCS density of states 
in S i at temperature Tj (i = 1,2), A^i(e,Ti) = 
A,(r-)/Ve 2 - A^T^efe 2 - A,(T,) 2 ], and e is the en- 
ergy measured from the condensate chemical potential. 
Furthermore, Aj(T^) is the temperature-dependent en- 
ergy gap, fi(Ti) = tanh(e/2fcsTi), 0(x) is the Heaviside 
step function, fee is the Boltzmann constant, and e is 
the electron charge. We note that both Qqp and Q^t 

vanish for T\ — T 2 [I3HT7] . while Q^t a ^ so vanishes 
when at least one of the superconductors is in the normal 
state, i.e., A^Ti) = 0. 

For a given the phases ip a and ipb are determined 
by the equation Ijs'm(p a = IjsixKfb which describes con- 
servation of the supercurrent circulating along the loop, 
where Ij^ oc R^Cb) is the Ambegaokar-Baratoff critical 
current [18] of junction J a (&). By defining r — Ij/Ij — 

Qint/Qlt = Q a JQ h qP = Rb/Ra (with < r < 1) as the 
degree of asymmetry of the SQUID junctions one gets 
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where x = $/$o- With the aid of Eqs. (4) and (5) Q int 
can be rewritten as 



Qint — Qint(Tl,T2)\ 
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which is analogous to the expression for the total Joseph- 
son critical current in a DC SQUID with generic junc- 
tions asymmetry |18j . In particular, for a symmetric 
SQUID (r = 1) we get 



Qint — 2Qi„ t (7l,T 2 ) 



7T$ 
$ 
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In the following we assume Ai(0) = A 2 (0) = A, with 
T c = (1.764fcs) _1 A representing the common critical 
temperature of the superconductors. 

Figure [2] (a) shows the interference component of the 
heat current Qint vs <& calculated for a few values of r at 
generic temperatures Ti and T 2 such that T 2 < Ti < T c . 
As it can be seen, Qi n t is a periodic function of $ which 
is maximized at integer values of <I>0: an d is modulated 
between the maximum given by Q int (l + r) and the min- 
imum given by Q| nt (l - r). For r = 1, Q in t(®) is thus 
modulated between 2Q b int and [see Eq. (|7j)]. By in- 
creasing r leads to a suppression of the modulation am- 
plitude combined with a reduction of the average value 
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FIG. 2. (Color online) (a) Interference heat current Qi nt 
vs $ calculated for a few values of r. Here we set generic 
temperatures Ti and T 2 such that T 2 < Ti < T c . (b) Total 
heat current Qtot vs <E> calculated for a few values of Ti at 
T2 = 0.1T C assuming r — 1. 



of the heat current. Eventually, the modulation ampli- 
tude is totally suppressed for r — 0, as only one junction 
is driving heat flow through the SQUID. Therefore high 
junctions symmetry is desired to maximize heat current 
modulation in the device. For that reason, in the follow- 
ing we will restrict our calculations to the case r = 1. 

Figure [2^b) shows the total heat current Qtot vs $ cal- 
culated at T 2 = 0.1T C for a few values of Ti. As expected, 
Qtot is $o"P er i oc iic and is minimized for integer values of 
$0 [note the sign minus in front of phase-dependent term 
in Eq. Q]. At the lowest T\ the total heat current is 
small, as Q qv and Qi n t are comparable in magnitude. By 
increasing Ti leads to the enhancement of both the aver- 
age heat current and modulation amplitude which origi- 
nate from larger temperature drop across the junctions. 
Further enhancement of Ti leads to a suppression of the 
heat current modulation amplitude which disappears at 
Ti = T c when Si is driven into the normal state. 

The above described behavior is emphasized in Fig. 3. 
In particular, the average Qtot over one flux quantuimdc- 

vs Ti for some values of T 2 . For each T 2 , (Qtot) mono- 

tonically increases by enhancing Ti, which stems from 
a larger absolute temperature bias across the SQUID. 

On the other hand, {Qtot^ becomes more and more sup- 
pressed at higher T 2 for fixed Ti. This suggests that 
a low T 2 combined with a substantial temperature bias 
are desirable to maximize the total heat current flowing 
through the interferometer. 

Figure [3] (b) shows the modulation amplitude 
5Qtot{T\)i defined as the difference between the maxi- 



fined as (Q to t) = %{Q qp - f Q\ nt ), is plotted in Fig. 
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FIG. 3. (Color online) (a) Average total heat current over 
one flux quantum ^Qtoty vs Ti calculated for some values of 

T%. (b) Total heat current modulation amplitude SQtot vs Ti 
calculated for the same Tb values as in panel (a). 



mum and the minimum values of Qtot > calculated for the 
same T 2 values as in panel (a). Specifically, for each T2, 
SQtot is a non-monotonic function of Ti, initially increas- 
ing and saturating at an intermediate temperature which 
depends on T 2 . In particular, by increasing T 2 , the max- 
imum of SQtot moves toward higher T\ values. Then, 
SQtot decreases at larger T\ vanishing at T c . The be- 
havior of Qtot therefore stems from the balance among 
several factors such as the temperature drop, the opera- 
tion temperature and the superconducting gaps. 

A practical experimental setup to observe these ef- 
fects is shown in Fig. [If a). The device can be realized 
with standard lithographic techniques, and consists of a 
tunnel-junction DC SQUID composed of identical super- 
conductors, and pierced by a magnetic flux <£>. Supercon- 
ducting leads tunnel-coupled to both SQUID electrodes 
and serving either as heaters (H) or thermometers (Th) 
allow to perturb and to accurately probe the quasiparti- 
cle temperature in the structure [I]. The superconduct- 
ing junctions provide nearly ideal thermal isolation of the 
SQUID electrodes [TJ and, therefore, we will neglect the 
thermal conductance through these probes. 

The relevant thermal model accounting for heat trans- 
port in the structure is shown in Fig. |4]jb) . Upon heating 
Si at temperature T\ the steady-state quasiparticle tem- 
perature T 2 will depend on the energy relaxation mecha- 
nisms occurring in S 2 . At low lattice temperature Tbath, 
typically below ~ 1 K, the predominant energy relax- 
ation mechanism in metals is related to electron-phonon 
interaction [TJ, Q e - p h(T 2 , Tb a th), which allows energy ex- 
change between quasiparticles and lattice phonons resid- 
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FIG. 4. (Color online) (a) Scheme of a practical setup con- 
sisting of a temperature-biased DC SQUID threaded by a 
magnetic flux Superconducting probes tunnel-coupled to 
both SQUID branches serve either as heaters (H) or ther- 
mometers (Th). (b) Sketch of the thermal model accounting 
for heat transport in the system. The arrows indicate the 
direction of heat currents for Tbath < T2 < Ti. 



ing at Tbath- Qe-ph in S 2 can be written as 
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coupling constant, and V is the volume of S 2 . At fixed 
Tbath, the steady-state T 2 (&) is obtained by solving the 
following thermal-balance equation 



Q tot (T 1 ,T 2 ,$) + Q e _ p , i (T 2 ,T i 



bath ) 



o. 



(9) 



In our thermal model we neglect phononic heat conduc- 
tion [T3], as well as heat exchange with the photonic 
environment due to poor matching impedance |211 122) . 
In the following calculations we assume an aluminum 
(Al) symmetric DC SQUID with bulk T c = 1.19 K, and 
R a = Rb = 1 kfi. As additional set of parameters we 
choose a S 2 electrode with volume V = 10 -19 m 3 , and 
E = 3 x 10 8 Wnr 3 K" 5 as typical for Al [TJ. 

The solution of Eq. (|9| is displayed in Fig. 5(a) where 
T 2 (<£>) is plotted for increasing Ti values at Tbath = 100 
mK. In particular, by increasing Ti leads to the enhance- 
ment of the average temperature T 2 over one flux quan- 
tum which stems from larger Qtot through the interferom- 
eter. Furthermore, while at low Ti the temperature mod- 
ulation amplitude <5T 2 (defined as the difference between 
the maximum and minimum values of T 2 ) is vanishing, 
at higher temperature it increases obtaining up to ~ 50 
mK at 800 mK. At larger Ti, 5T 2 tends to decrease due 
to temperature-induced suppression of the energy gaps 
in both superconductors. We note that this behavior re- 
flects basically the same picture provided by Qtot [see Fig. 
[5Jb)]. It is worthwhile to mention that by reducing the 
S 2 volume so to suppress Q e - P h does not yield necessarily 
enhanced <5T 2 for fixed T\ and Tbath, as in such a case the 
temperature drop across the SQUID would decrease. In 
the opposite situation of larger V, enhanced Q e -ph would 
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FIG. 5. (Color online) (a) Quasiparticle temperature T2 vs 
$ calculated for a few values of 21 at Ti, a th = 100 mK. (b) 
Flux-to-temperature transfer function T vs $ calculated for 
selected values of Ti at T bat h = 100 mK. (c) T2(<1?) calculated 
for different T bath at Ti = 1 K. (d) T($) calculated for the 
same values as in panel (c). 



suppress <5T 2 as well. Analogously, by largely increasing 
R a leads to a reduced Q to t and to a suppression of 6T2, 
whereas for low R a reduced temperature drop across the 
junctions would prevent large 8T2 values. 

In analogy to a conventional electric SQUID, we can 



define for the heat-flux modulator a figure of merit in 
the form of a flux-to-temperature transfer coefficient, 
T(<S>) = dT 2 /d<P, which is plotted in Fig. 5(b) for se- 
lected values of Ti at T bath = 100 mK. Specifically, T(<J>) 
behaves similarly to T2($) as a function of Ti, and ob- 
tains sizable values as large as ~ 125mK/$o a t 800 mK. 

The role of bath temperature is shown in Fig. 5(c) 
which displays T2($) calculated for increasing Ttath val- 
ues &t Ti = 1 K. An increase in T^/j leads to a smearing 
of T2($), and to a reduction of ST 2 due to the combined 
effect of enhanced Q e - P h, reduced temperature drop and 
temperature- induced gap suppression in S2. As expected, 
full suppression of 5T 2 occurs at T^y, = Ti, as in such 
a case Si and S2 are at the same temperature. This be- 
havior directly reflects on T(<1>) which is shown in Fig. 
HJd) for the same values as in panel (c). 

We shall finally comment the temperature modulation 
speed which could be achieved in the device. The latter 
is mainly limited by the relaxation time r required by 
quasiparticles in S2 to thcrmalize with lattice phonons, 
since the R a C time constant of the SQUID junctions (C 
is the capacitance) can be made much smaller than r 
with a suitable choice of parameters. In particular, in 
the 0.5 ...IK temperature range, r -1 is of the order of 
~ 1 ... 10 MHz for Al [53], whereas at lower Tb at h it is 
drastically reduced owing to increased electron-phonon 
relaxation time [33] 121] • We stress that the modulation 
frequency could be increased by more than a factor of ten 
by exploiting superconductors with enhanced electron- 
phonon coupling like, for instance, tantalum [24] . 

In summary, we have proposed and analyzed a phase- 
controlled superconducting heat-flux quantum modula- 
tor. In an Al-based design, it can provide temperature 
modulation amplitudes up to 50 mK and transfer func- 
tions as large as ~ 125 mK/$ below 1 K. This struc- 
ture appears potentially attractive for the design of novel- 
concept thermal devices operating at low temperature. 
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